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CRACK PROPAGATION AND ARREST . 

IN PRESSURIZED CONTAINERS ^ ^ 

by 

F. Erdogan, F. Dolale, a»»d J. A. Owc^arek 
Lehigh University, Bethlehem, Pennsylvania 

AbiPuxct . The problem of crack propagation and arrest in a finite volume 
cylindrical container filled with pressurized gas is considered . It is 
assumed that the cylinder contains a symmetrically located longitudinal 
part-through crack with a relatively small net ligament. The net ligament 
suddenly ruptures initiating the process of fracture propagation and de- 
pressurization in the cylinder. Thus che problem is a coupled gas dynamics 
and solid mechanics problem the exact formulation of which does not seem 
to be possible. It is formulated by making two major assumptions, namely 
that the shell problem is quasi-static and the gas dynamics problem is 
one-dimensional. The problem is reduced to a proper initial value problem 
by introducing a dynamic fracture criterion which relates the crack accel- 
eration to the difference between a load factor and the corresponding 
strength parameter. The main results are demonstrated by considering two 
examples, an alvuninvtm cylinder which may behave in a quasi-brittle manner, 
and a steel cylinder which undergoes ductile fracture. The results in- 
dicate that generally in gas-filled cylinders fracture arrest is not pos- 
sible unless the material behaves in a ductile manner and the container is 
relatively long. 


***^This work was suppc rted by the Materials Division, NASA-Langley under 
the Grant NGR 39-007-011 and by NSF under the Grant ENG 73-045053 AOl. 
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INTRODUCTION 


In studying the fracture failure of relatively thin-walled pressure 
vessels and pipes, in general one may consider the phenomenon in three 
different phases. The first is the subcritical crack propagation starting 
from a localized imperfection. This initial flaw which may form the nu- 
cleus of a dominant crack in the cylinder wall may be a material or manu- 
facturing defect, or may be accidentally induced. Even though such flaws 
are just as likely to be fully imbedded interior defects, from the view- 
point of subcritical crack growth, it is the surface flaw which is of more 
critical concern. The reason for this is the presence of higher level 
stresses near and at the surfaces and the likelihood of ready exposure to 
adverse environmental conditions. Once the crack is nucleated, the sub- 
critical crack propagation may be due to low cycle fatigue and stress 
corrosion cracking. 

In the heavy-section low toughness structures if the subcritically 
growing crack reaches a "critical" size, it is possible to have a plane 
strain type of fracture leading to catastrophic failure. However, in 
most cases, particularly in relatively thin-walled cylinders, there is 
sufficient amount of plastic deformation around the leading edge of the 
part-through crack so that before the crack reaches the critical size 
the net ligament between the crack front and the opposing surface would h 
fully yielded (Figure 1). In such cases, the second phase of the failure 
process would then be {generally ductile) fracture of the net ligament. 

The net ligament rupture may occur as a result of progressive stable crack 
growth. However, for materials with relatively high fracture toughness, 
perhaps a more likely rupture mechanism is the "net ligament plastic in- 
stability" described in [l]. 

Even though in some cases it is theoretically possible to have "leak" 
before "burst", in practical situations it turns out that in both brittle 
as well as ductile type of fracture, the internal pressure which is high 
enough to cause the net ligament to rupture is usually above the critical 
level which can be sustained by a cylinder with a through crack having 
the same length 2a^ as the part- through crack [l] (Figure 1) . Thus, the 
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1 ,iial phase of the fracture failure in pressurized cylinders is the prop- 
agation of a through crack. This is clearly a coupled fracture dynamics 
and fluid mechanics or gas dynamics problem. Wiether the fracture propa- 
gation would result in catastrophic failure or crack arrest would, of 
course, depend on the nature of decompression in the cylinder as the 
liquid or gas is discharged through the orifice provided by the propagating 
crack . 

From the structural view point the tools dealing with the subcritical 
crack growth phenomenon appear to be fairly well-understcod and quite ad- 
equate [ 2 , 3 ]. If the particular part-through crack problem can be solved 
(see, for example [4,5]), then the crack propagation rate under fatigue 
or stress corrosion cracking may be estimated by using certain empirically 
established models [2, 3, 6, 7]. A quantitative modeling of fully-yielded 
net ligament in shells is somewhat more complicated. However, in this 
case too adequate results may be obtained by properly considering the 
plasticity effects [l,8,9]. On the other hand, it would be no exaggera- 
tion to state that the crack propagation phase of the fracture failure 
following the net ligament rupture is not understood at all. In this paper, 
after making some remarks regarding the fracture dynamics in relatively 
thin-walled pressurized containers mostly for the purpose of pointing out 
the extremely complex nature of the problem, the quasi static fracture 
propagation problem for the finite volume containers is considered. Two 
types of fracture models applicable to materials with essentially tie 

or ductile fracture behavior are developed and some numerical exampiet:- 
are given. 

2. ON THE FRACTURE DYNAMICS IN PRESSURIZED CONTAINERS 

Consider a relatively thin-walled pressurized cylindrical container 
of finite volxune. Let a symmetrically located meridional through crack 
of length 2a^ suddenly appear in its wall as a result of net ligament rup- 
ture (Figure 1) . The ensuing phenomenon may be considered as an initial 
value problem with, for example, the components of the displacement vector 
Uj^, u^, and u^ in the cylinder wall and the components of the velocity 
vector Vj^, v^, and v^, the pressure p, the temperature T, and the density p 
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of the liquid or gas inside the cylinder as the basic state variables 
which are functions of time and space coordinates. Considering the sym- 
metry of the problem and assigning that no crack branching takes place, the 
problem has at least one more unknown function which is the crack length 
a(t) . The crack surface displacements which determine the size and the 
shape of the "orifice" are given by the displacement components Uj^, u^, 
and Uj. By using the basic physical principles, theoretically it should 
be possible to write nine field equations (e.g., equations of state and 
energy, and conservation of mass and momentum) , to account for the nine 
basic unknown functions u^, v^, (i=l,2,3), p, T, and p. These equations 
may be expressed in terms of a system of differential operators as follows: 

L^(u^,v^,p,p,T) =0, (i=l,..,9, j=l,2,3) . (1) 

Equations (1) will be subject to standard initial conditions, the condi- 
tions of continuity on the inner surface of the cylinder regarding the 
traction and the velocity components, and the traction boundary conditions 
on the outer surface of the cylinder and on the surface of the crack. Even 
if there is no crack propagation after the rupture of net ligament, the 
general dynamic problem is intractable. The reason for this is that in 
considering the gas dynamics or the fluid mechanics, it does not seem to 
be possible to express the boundary conditions at the opening on the cyl- 
inder wall created by the crack in an exact manner. 

If the crack continues to propagate after the rupture of net ligament, 
the problem is further complicated by the fact that it now has a moving 
boundary. In this case, technically one may complete the formulation of 
the problem by considering an additional unknown function which represents 
the crack size, such as the half crack length a(t) , and an additional con- 
dition, such as a fracture criterion of the form 

F(u.,p,a,Cj) > 0 (2) 

which states that a certain minimum condition must be met for propagation 
of the crack. In (2) the displacements u^, (i=l,2,3) represent the stress 
or the strain field in the cylinder, the nonuniform pressure p represents 
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the applied load which may appear in the criterion explicitly as well as 
through u^, and (j=l,2,,,) are the constants representing the yield 
and fracture resistance behavior of the material. Ordinarily, in struc- 
tural mechanics the form of the fractiire criterion as given by (2) and, 
for example, expressed in terms of the stress intensity factor k > kj^, 
the strain energy release rate G > G^^, the J-integral, or the Kj^-curve 
is adequate if one is only interested in the question of stability of the 
crack under a given set of loading conditions which are usually either 
fixed load or fixed displacement type. However, it is clear that in a 
coupled problem of fracture dynamics-fluid mechanics, it is not sufficient 
to say that the fracture will take place if certain minimum conditions are 
met and exceeded. It is at the same time necessary to make some specific 
statements regarding the dependence of the fracture acceleration and de- 
celeration or fracture velocity on the load level. More specifically, 
the fracture criterion must express the dependence of some time rate of 
change of the fracture area on the difference between a load factor repre- 
senting the intensity of the applied load arid a corresponding strength 
parameter which represents the fracture resistance of the material . Lack 
of phenomenological understanding and adequate modeling has been one of 
the major difficulties in studying the problem of fracture dynamics in 
structural solids. Because of this in most of the existing studies rather 
than considering a physically acceptable fracture criterion, an inverse 
technique is adopted by specifying the crack velocity [5,lo]. 


3. A QUASI STATIC MODEL FOR GAS-FILLED CYLINDERS 

Consider now a cylindrical shell of length Z, mean radius of curvature 

R , and thickness h. Initially let the cylinder be filled with gas of 
o 

pressure p^ and {absolute! temperature T^. It is assumed that the cylinder 
contains a symmetrically located part-through longitudinal crack with a 
relatively small net ligament thickness {Figure 1) . At time t=0 the net 
ligament ruptures creating a through crack of length Let the pressure 

p be greater than the value which can be sustained by a cylindrical shell 
having a through crack of length 2a^. Thus, the crack will start growing 
upon the rupture of net ligament and, because of the gas outflow through 
the crack opening, the container will start to depressurize. From the 
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view point of crack propayation, sitnilar to flat plates, under constant 
load the cracked cylindrical shell is basically unstable. However, be- 
cause of depressurization there is a pos ibility of crack arrest before 
catastrophic failure. It is therefore clear that the problem of fracture 
dynamics is strongly coupled with the problem of gas dynamics. 

To develop any type of an acceptable analytical model for the problem 
two major obstacles need to be overcome. As pointed out earlier, the 
exact formulation of the related gas dynamics problem does not seem to be 
possible. Therefore, even if the necessary mathematical techniques to 
solve such problems were to be available, without any approximating assump- 
tions the problem would be intractable. This then is the first major as- 
pect of the phenomenon which requires certain physically acceptable simpli- 
fying assumptions so that an analytical treatment of the problem is 
possible. The second muior aspect of the phenomenon refers to fracture 
dynamics which is fundfimentally much less understood. In the present 
coupled problem to use the so-called inverse method which ass’umes the 
velocity or the acceleration of the crack beforehand would be meaningless. 
The key question in this respect is this; in v.»hat way is the crack prop- 
agation behavior (say, its velocity or acceleration) dependent on the load 
factor which represents the intensity of the applied load and the geometry? 
Without making some statement to this effect, clearly the crack propagation 
and arrest problem cannot be treated as a dynamic pVienomenon. 

Consider now the first aspect of the jjhenomenon, namely the formula- 
tion of the coupled solid mechanics-gas dynamics problem. To render the 
problem tractable and bring it down to manageable proportions, in this study 
the following approximating assumptions will be made; 

(a) In the shell analysis the inertia effects will be neglected and 
the problem will be treated as being quasi-static. Even though there are 
no available solutions to the dynamic shell problem, using the flat plate 
results as guide it may be stated that the inertia effects will be neglig- 
ible provided the crack velocity remains below approximately one-fourth of 
the shear wave velocity of the shell material [lo]. As will be seen from 
numerical examples, in the finite volume containers this ceiling for the 
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crack velocity and hence# the quasi-static assumption does not seem to be 
very unrealistic. 


(b) The compressed gas in the contaiier behaves like an ideal gas 
with constant specific heats. 

(c) The pressure drop in the container is sufficiently slow so that 
the gas dynamics problem may bo treated as being "one-dimensional*' with 
the (space-independent) basic variables pressure p# temperature T, and 
density p. 


(d) The discharge through the crack opening is at sonic velocity with 
a constant discharge coefficient Cj^, i.e., it is assumed that 


P 

*^e 

P(t) 




(3) 


where p is the gas pressure in the containers, p^ is the environmental pres- 
sure, Y“Cp/c^# Cp and c^ being the specific heats under constant pressure 
and constant volume, respectively. 


(e) The crack opening or the "orifice" may be approximated by an 
ellipse with the crack length and the crack opening displacement in the 
center of ihe crack as the major axes. Using the results given in [n] 
and extrapolating them the crack opening area may be approximated by 

2R , .^2 

D,t, = 


(4) 


where E is the elastic modulus of the cylinder. 


Under these assumptions it is seen that the unknown functions are p, 
p, T, and a with t as the independent variable. The conditions to deter- 
mine these functions are 

(i) the equation of state for the gas 

^ = RT (5) 

P 

where R is the gas constant; 

(ii) the isentropic process equation 

— Y 

pp = constant; (6) 
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(iii) the conservation of mass which states that the rate of mass 
outflow through the orifice is the same as the rate of mass reduction in 
the container and which may be expressed as 


3t ** • 


where is the initial mass of the gas, V is the container volume, is 
the discharge coefficient, and D is the crack or orifice area. The quanti- 
ties 0 ^ and u* are, respectively, the gas density and velocity at the small- 
est exit section where sonic flow conditions are assumed to prevail, and 
are given by 

1 

0. - - lily "I')'’ • 18) 


Eliminating p and T, from equatio.iS (4) to (8) it follvjws that 




(Y-D/2Y hE 


1.47a’. 2^(5 y-1)/2y 


where T and p are the initial gas temperature and pressure in the container, 
o o 

(iv) Dynamic fracture criterion. A fracture criterion stated in the 
form of an inequality such as (2) can only refer to fracture initiation or 
to a static condition. In order to treat dynamic crack propagation prob- 
lems, the criterion must also be dynamic in nature and roust properly relate 
"the crack driving force" to the crack length and its time derivatives. In 
expressing the dynamic fracture criteria used in this study, two basic 
physical assumptions are made which are intuitively motivated. It is as- 
sumed that first the crack acceleration or velocity is a function of the 
"crack driving force" and secondly, if the crack driving force is negative, 
in structural metallic materials there is a fracture deceleration phase and 
the deceleration rate is much higher than the crack acceleration rate. The 
next ?tep is the proper selection of the "crack driving force" which will 
be dependent on the fracture behavior of the material. In brittle or quasi- 
brittle type materials it is realistic to assume that k-k^ or G-G^, is 

TST) . 

Here the critical values k(. and may be velocity dependent and may not 

be the same as the plane strain values Kjq and Gj^. In this sense, the model 
has some flexibility and may accommodate the phase of the fracture which may 
not be completely plane strain in nature. 
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the crack driving force where k is the time -dependent stress intensity 
factor at the crack tip, k i& a critical value which has to be exceeded 
for fracture initiation and gr-'*wth, G and are the strain energy re- 
lease rate and fracture toughn' ss and are related to k^^ k^, respectively. 

The experiments on glass plates indicate that "the crack has no in- 
ertia". The meaning of this statement and the related conclusions are 
that the crack velocity is always in phase with the crack driving force, 
the velocity is a functif^n of the crack driving force alone, and crack 
branching occurs when the input energy (or the crack driving force) becomes 
too high to be dissipated by a single crack and takes place at a certain 
limiting crack velocity (approximately half the shear wave velocity) . 

Since the limit iny crack velocity implies increased fracture resistance 
at higher crack velocities, for ideally brittle materials the dynamic 


fracture criterion may be expressed by relating the crack velocity to 


k-kj, or G“G^. Assuming that G'^k^, one may, for example, write 

[o , for k < k^ 


(10) 


where dj^,..,dj^ are constant, k is generally a function of a, and k^ is a 
function of da/dt. In the range of crack velocities for which the inertia 
effects are negligible, k^ may be assumed constant and the differential 
equation (10) with (9) would determine the functions a(t) and p(t) . 

Intuitively it could be argued that in metallic materials with varying 
degrees of plastic behavior the instantaneous change in velocity implied 
by (10) is not possible and the phenomenon must exhibit some inertia. In 
this case the dynamic fracture criterion must include crack acceleration 
as well as the velocity. Also, if the plasticity effects are not apprec- 
iable, k or G may be used as the crack driving force and the fracture cri- 
terion in its simplest form may be expressed as 


Model I; 

^2^ [bi(k-kc) , for k > k^, 
"\bj^(k-k^) - b2(k'k^)^ 


for k < k , 
c 


( 11 ) 
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or Modal II: 

A ^1*= ->'cl 

|b (k2-k^) 

' 4 C 

where b. are constant. 

1 4 

The stress intensity factor k 
of the shall problem (e.g., [l.?,13]). Extrapolated for large values of 
the shell parameter X, the stress intensity factor for the pressurized 
cylindrical shell with an axial crack may be approximated by 


k “ v^TTa A(X) , 

(13) 

fo.5X +0.6 + 0.4e~^*^^ for X < 5 


A(X) e< 

(14) 

k.761 A -1.9 , for X>5 , 


X = [12(1-V^) j"* a//Rji 

(15) 


where v is ♦-f., Poisson's ratio. 

In the case of ductile fracture one may need a fracture criterion 
which takes into account the ductile behavior of the material. Two such 
criteria based on the crack opening stretch and plastic instability were 
described in [s] and [l]. in either criteria the material’s fracture re- 
sistance parameter may be expressed in the form of a "resistance pressure" 
wfich is a ftinction of the crack length. For example, for steel cylinders 
the resistance pressure pj^ or the load carrying capacity of the vessel may 
be expressed as [l] 

{l-n)c ] (16) 

where n and a are constant (0<n<l, a>0) , X is the shell parameter given 
by (15), and 

ha 

Py = ~ , Oy = OygU+6) , (0.05 <B< 0.15} . (17) 

o 

Here is the yield strength of the material, Oy is defined as the flow 
stress which is somewhat higher than O , and p is the pressure correspond- 

Xo y 

ing to a "fully-yielded" cylinder. Thus, in this case the "crack driving 

. 2 2 

force may be taken as either P~Pj^ or p -p^ and the following dynamic fracture 


, for k > kg , 

- ^ . for k < k^ 


( 12 ) 


is obtained from the elastic solution 
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vrriterla may be expressed: 


Model III: 

, for p>pj^ , 

Ui(P“Pr> ~ c^Cp-p^)^ , for p<pj^ , 


or Model IV; 

^ ^ ^3(p^-Pr) . for p>pj, . 

dt2 % 3 (p 2 -p 2 ) . c^(p% 2)2 ^ f or p < p^ , 


( 18 ) 


(19) 


where again are constant. These as well as the constants 

have to be determined experimentally. Such experiments may be conducted 
on flat plates under ideal leading conditions and crack geometry. 

The coupled depressurization and quasi-static crack propagation prob- 
lem may now be solved by considering (9) and any one of the equations (11) , 
(12) , (18) or (19) . The problem is a typical initial value problem which 
is highly nonlinear. It may be solved n’omerically by defining 

~ = v(t) (20) 

and by observing that the differential equation (9) and the Models I to IV, 
i,e,, (11), (12), (IS) and (19) are of the form 


f(a,p) , 


d a dv , , , . , . . 

— r = TT =* 9, (a»P) * (1=1,., ,4) 

dt^ dt 1 


Using the initial conditions 


a(0) = a , v(0) = 0 , p(0) = p^ 

o o 


one may then write 


p(t) = Pq + 


a(t) = a + 
o 


v(t) 


= i: 


f(a,p)dt , 
v(t)dt , 

p)dt , (i=l,..,4) 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 
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and use any one of the integration fumuiab to obtain the unknown functions 
y, a, and v. The numerical calculations are terminated when either 

v(t^) - 0 , a(t.) <0.el/2 , (27) 

or 

a(tp) - 0.8^/2 , v(tp)>0 . (28) 

v'hev-e S. is the total length of tho cylindrical container. The former situa- 
tion corresponds to fracture arrest t^ bning the arrest time and the latter 
basically implies catastrophic failure. The approximate limit a«0.8V2 is 
imposed by the shell analysis which is adopted from the infinitely long 
cylinder solution. 

4. EXAMPLES AND DISCUSSION 

P.s an example for a pressure tank with brittle or quasi-brittle frac- 
ture behavior a 20 in. diameter aluminum cylinder and for a pressure tank 
with ductile fracture behavior a J6 in. diameter steel cylinder is consid- 
ered. In both cases the gas is assumed to be a biatomic ideal gas. The 
relevant information regarding the numerical data used in the calculations 
is given in Table I. The initial value problem as formulated by equations 
(9) and (10) or (20) to (23) was solved by using two sets of values for 
the crack propagation constants b^, c^, (i>l,..4) and d^ (i.e., N«1 in (10)), 
which are shown in i''\ble II. Other variables in the solution were the 
length of the cylinder Z, the discharge coefficient C^, and the amount of 

the initial overload such as k(0)-k or p -p (0) . Some of the typical re- 

O O R 

suits are shown in Figures 2 to 15. 

Figures 2 to 5 show the results for the aluminum cylinder which was 

assumed to behave in a quasi-brittle manner. Here the length was varied 

between 20 in. and 200 in., the discharge coefficient was 0.5 < < 0.8, 

and the initial pressure overload was such that 50 < k(0)-k^ < 500 psi /In. 

Models I and II were used to calculate a(t), p(c) , v(t) , and k(t) using the 

constants shown in Table II. Figures 2 and 3 show the results obtained 

from Model I, i.e., equation (11), by using the constants given in columns 

A and B of Table II, respectively. In these as well as in the subsequent 

figures the computed values p, a, v, k, and p are normalized with respect 

R 
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Table I. Data used in the examples 



Aluminum (Models I, II) 

Steel (Models III, IV) 

a 

0.5 in . 

0.0127 m 

1.9 in. 

0.04926 m 

o 





h 

0.1 in. 

0.00254 m 

0.403 in. 

0.01024 m 

R 

10 in. 

0.254 m 

18 in. 

0.4572 ra 

o 


5.08 m 

r 50 in. 
^500 in. 

r 1.27 m 
^12.70 m 


200 in. 

E 

10^ psi 

10 2 
6.895 X 10 N/m‘ 

3 X 10^ psi 

20.685 X 10^°N/m^ 

V 

0.3 


0.3 



70 ksi 

48.265 X lO^N/m^ 

65 ksi 

44.8175 X lO^N/m^ 

0 



75 ksi 

51.7125 X lO^N/m^ 

Y 





k(0) 

45.05 ksit'Hi 

7 -3/2 

4.945x10 Nm 



k 

45 ksi/In 

6.94 X lO^Nm”^'^^ 



c 





Po 

' -i.65 psi 

2.1 X loS/m^ 

1108.53 psi 

7643.314 X lO^N/m^ 

PrfO) 



p^-10 psi 

7574.364 x lO^N/m^ 

R 

238.752 


238.752 


T 

500®R 

278°K 

525®R 

292®K 

o 





Y 

1.4 


1.4 


S 

,0.5 

^0.8 


rO.5 

'o.8 






Table II. Crack propagation constants (iji 
units of lb, in, s<» ) 


Model 

■ 

A 

D 

I 

H 

8 X lo"^ 

8 X lO"^ 

Eq. (11) 

^2 

4 X lo"*^ 

4 X lO"® 


H 

8 X 10“^° 

-9 

8 X 10 

II 

Eq. (12) 

H 

4 X lo"^"^ 

4 X 10 



0.002 

0.02 

III 




Eq. (18) 

^2 

0.001 

0.01 



0.4/(p^-p^(0)) 

4/(p^-p^(0)) 

IV 

Eq. (19) 

^4 

80/(p^-p^(0)) 

800/(p^-p^(0)) 

Eq. (10) 

^1 

O 

1 

CO 

10-'^ 


to a set of constants p , a , v , k , and p„ which are different for each 

n n n n Rn 

example and are given by the figure captions. Thus the quantities shown in 
the figures are defined by 

a = a(t)/a^ , p = p{t)/p^ , v = v(t)/v^ 

k = k{t)/k^ , PR=PR(t)/Pj^ . (29) 

From the crack arrest view point these figures correspond to the most favor- 
able conditions, that is, they were calculated for longer length (£= 200 in.), 
greater discharge coefficient (Cj^ = 0.8), and smaller overload (k(0)-k^=50 
psi v^in) . Even in this case, in spite of the sharp decrease in the pressure 
and the stress intensity factor after a certain time, at a =0.8 £/2 the 
crack velocity v was still increasing. The basic conclusion here is that in 
gas-filled cylinders behaving in quasi-brittle manner it is unrealistic to 
expect any crack arrest. 
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Figures 4 and 5 show the results obtained from Model II * i.e., equa- 
tion (11), by using the crack propagation constants given in columns A and 
B of Table II. Note th qualitatively these results are quite similar to 
those obtained from Model I. 

The results for the steel cylinders are shown in Figures 6-13. Fig- 
ures 6-9 show the results obtained from Model III, i.e., equation (18) for 
various combinations of Z, and c^^ and c^ (see columns A and B in Table 
II) . Here the figures show the driving pressure p and the resistance 
pressure p_ with the same scale. It is seen that for the short cylinder 
(■C=50 in) before the pressure drops below p , i.e., before t)ie crack starts 
to decelerate, the crack length reaches 80 percent of the cylinder length, 
hence no crack arrest appears possible (Figures 6,7) . On the other hand, 
in the long cylinder (£= 500 in) the crack is arrested before the crack 
length reaches £/2. Changing various constants in the calculations indicated 
that in the ductile fracture model, from the crack arrest view point the 
most important variable was tlie cylinder length. For example, comparing 
Figures 8 and 9 it may be seen that in spite of one order of magnitude dif- 
ference in crack propagation constants c^, c^ (see columns A and B in Table 
II), the results are qualitatively the same. Higher values of c^^ and c^ 
simply shortens the crack arrest time. In the results shown in Figures 8, 

9, 12, and 13 where crack arrest took place, the calculations were stopped 
when the crack velocity dropped to zero. 

The results obtained form Model IV (equation (19)) are shown in Fig- 
ures 10-13. Except for the shortened time scale, these results are quite 
similar to those obtained from Model III. 

The results obtained from the velocity model (10) with N = 1 are shown 
in Figures 14 and 15. In these calculations too the influence of crack ve- 
locity on k^ or is neglected. The highest velocity observed here is 
approximately one- fourth of the shear wave velocity of the material, meaning 
that the inertia effects may still not be all that crucial. It is seen that 
qualitatively these results are quite similar to those obtained by using the 
acceleration Models I and II, the important difference being in the time 
scale. In t>oth cases there does not seem to be sufficient time for any 
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significant depressurisation. Hence, the probl<an is analogous to the cen- 
trally cracked flat plate problem under constant tensile stress with the 
additional magnification in the "crack driving force" k coming from the 
shell curvature effect. 

The saR^le calculations given in this paper indicate that in order to 
deal with the crack propagation and arrest phenomenon the fracture cri- 
terion must be cast in dynamic form relating the time rate." of the crack 
length to the overload. Developing such dynamic fracture criteria does 
not seem to be a major problem. The main difficulty in this respect lies 
in determining the related crack propagation constants. This may be done 
by conducting carefully designed experiments in %<hich the crack length and 
the appropriate load factor are accurately measured or calculated and the 
corresponding resistance parameter is independ ntly estimated. The pre- 
liminary calculations given in this paper also indicate that in most finite 
volume containers it is very likely that the crack velocity will always be 
'ufficiently low to justify the quasi-static assumption in the shell 
analysis. 
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Figure 1. Geometry of the pressurized cylindrical container. 
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Figure 2. Results fui the 20 in. diameter aluminum vessel; £- 200 in 
Cq- 0.0, Model I with constants in column A, Table II. Th< 
normalization constants: 300 psi, a„ ■ 80 in., * 80 

in/sec., ICj, ■ 16 X 10^ psi /In. 
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Figure 3. Results for the 20 in. diameter aluminum vessel; £- 200 in 
C[)-0.8, Model I with constants in column B, Table II. Thi 
normalization constants: p,^ - 300 psi, a^ - 80 in., - 300 

in/ sec. kj) - 2 X 10® psi /IrT. 





Figure 4. Results for the 20 in. diameter aluminum vessel; £ ■ 200 in., 
Cd- 0.8, Model II with constants in column A, Table II. The 
normalization constants: ^* 300 psi, a^ • 60 in., ■ 400 

in/sec., 20x10^ psi 
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Figxire 6. Results for the 36 in. diameter steel vessel; £*50 in., 

Cjj“0.5, Model III with constants in column A, Table II. The 
normalization factors: ■ 1000 psi, a,^ ■ 20 in., ■ 7 in/sec. 

PRn * 1000 Psl* 





Figure 7. Results for the 36 in. diauneter steel vessel; t*50 in., Cjj*0.5, 
Model III with constants in column B, Table II. The normalization 
factors: Pn “ * 1000 psi, ■ 20 in., “ 20 in/sec. 
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Figure 8. 
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Results for the 36 in. diameter steel vessel; £* 500 in., 
Cq<- 0.8, Model III with constants in column A, Table II. 
Normalization factors: ° pj^yj » 1000 psi, a„ ■ 65 in., 10 

in/sec. 
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Results for the 36 in. diameter steel vessel: £= 500 in., 
Cn=0.8, Model III with constants in column B, Table II. Nor- 


malization factors: 


1000 psi, a » 75 in., ■ 36 in/sec 





Figure 10. Results for the 36 in. diameter steel vessel; £ « 50 in., Cp-0.5, 
Model IV with constants in column A, Table II. The normalization 
factors: ■ 1000 psi, ajj * 20 in., ■ 25 in/sec. 
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Figure 12. Results for the 36 in. diameter steel '.essel; £* 500 in., 

Cq* 0.8, Model IV with constants in column A, Table II. The 
normalization factors; ■ 1000 psi, a^^ » 30 in., 

V » 40 in/sec. 
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Results for the 36 in. diameter steel vessel; £* 500 in., 
Cjj*0.8, Model IV with constants in column B, Table II. 

The normalization factors; * 1000 psi, aj^ » 100 in., 

v_ * 150 in/sec. 
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Figure 14. Results for the 20 in. dicuneter aluminum cylinder; £* 200 in., 
Cjj»0.8. The velocity model, Eq. (10), d. * 10"®; normalization 
factors: a_ * 50 in., p * 300 psi., v„* 30,000 in/sec., 

k„-2xio® psi /nr. 




